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Article Info ABSTRACT 
Article history: The genetic algorithm (GA) is an adaptive metaheuristic search method based on the 
Received Aug 3, 2021 process of evolution and natural selection theory. It is an efficient algorithm used for 


solving the combinatorial optimization problems, e.g., travel salesman problem (TSP), 
linear ordering problem (LOP), and job-shop scheduling problem (JSP). The simple 
GA applied takes a long time to reach the optimal solution, the configuration of the GA 
parameters is vital for a successful GA search and convergence to optimal solutions, it 
Keywords: includes population size, crossover operator, and mutation operator rates. Also, very 
recently, many research papers involved the GA in coding theory, In particular, in 
the decoding linear block codes case, which has heavily contributed to reducing the 
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Automorphism group 


Crossover : complexity, and guaranting the convergence of searching in fewer iterations. In this 
Genetic algorithm paper, an efficient method based on the genetic algorithm is proposed, and it is used 
Low density parity check for computing the Automorphisms groups of low density parity check (LDPC) codes, 
codes the results of the aforementioned method show a significant efficiency in finding an 
Mutation important set of Automorphisms set of LDPC codes. 
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1. INTRODUCTION AND PRILIMINARIES 

There are varying methods in coding theory which addresses its application, one of them is through de- 
termining the Automorphisms groups of codes, they allow us to determine the structure of the codes, classifying 
them and help the decoding algorithm. This remains a challenge since determining the whole automorphisms 
groups of codes is difficult, except finite simple groups which have been realized using the sporadic groups 
(e.g, the aumorphism group of golay codes are mathieu groups). 

Recalling that the hamming distance between any two codewords (vectors) c, c’ in F is defined to be 
the number of coordinates in which c and c’ differ. A binary linear [n,k,d]-code C over Fə is a k-dementional 
subspace of the vector space F3, where: 

: ' P 
d=d(C) = ae d(c,¢c) = re wt(c) (1) 
and its generator matrix G is a k x n matrix whose rows is the basis of C. 

Let C be a binary linear code and G its generator matrix, considering the action of the symmetric 

group Sn on the G columns. For all ø in Sn, denote by G.o the matrix obtained from the permutation of the G 
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columns. Let o € Sn, € = (c1, C2, ..., Cn) E C: 
a(c) = o (c1, C2,- » Cn) = (Co(1): Co(2) +- -3 Co(n)) (2) 
a(C) = {a(c),c E€ C} (3) 


any permutation of the G columns which maps the rows of G into rows of the same matrix, is called an 
automorphism of C. The set of all automorphism permutations forms a subgroup of Sn, denoted by Aut(C'): 


Aut(C) = {0 € C,a(C) = C} (4) 


let A be a group, A is an Automorphism group of C if A C Aut(C) and A is the Automorphism group of C if 
A = Aut(C) B]. 

This paper, mainly focuses on the computation of the automrphisms groups of LDPC codes. The sec- 
tion 2, includes some definitions, details, also it presents related works using genetic algorithm (GA). In section 
3, the GA-based method is described, including the fitness function, stochastic crossbreeding, and stochastic 
operators. The results are presented in section 4. Section 5 is devoted to the conclusion and perspectives. 


2. RELATED WORKS 
2.1. Low density parity check codes 

Gallager devised the low density parity check (LDPC) codes, often known as Gallager codes, in 1962, 
they are class of linear block codes, defined by sparse parity check matrices, where each column contains a 
small fixed number w, of ls and each row contains a small fixed number w, > we of ls [B]. Due to the 
limited characteristics of computers at that times, this class of linear code was absent till 1990s where they 
have been reinvented through the Macky and Neal works, its has been shown that LDPC codes performance 
is near to Shannon limit performance with belief propagation algorithm (BPA) [4]. There are characteristics 
that distinguish LDPC codes from Turbo codes, such as superior performance when the block length is large, 
enormous flexibility, easy description and subsequent theoretical venerability, decreased decoding complexity, 
and so on [5]. 

There is an algebraic representation, the LDPC code is denoted as (n, we, wr), where n is the binary 
linear code length, we is the number of 1s in the column of the sparse parity check matrix (i.e. the column 
weight), and w, is the number of 1s in the row in of the sparse parity check matrix (i.e. the row weight) as 
illustrated in Figure 1, if w. and w, are invariant, it’s called regular LDPC codes, else it’s called irregular LDPC 
codes. Both of the two must satisfy this following condition: 


cH’ =0 (5) 


where c is a codeword and H is the sparse parity check matrix. There is another representation for LDPC codes 
which is trough Tanner graphs (graphical representation of the sparse parity check matrix), they contain two 
class of nodes, variables nodes, they represent the sparse parity check matrix columns, and check nodes, they 
represent the sparse parity check matrix rows. for each nonzero hij of H, an edge will be presented between 
check node i and variable node j as illustrated in Figure 2. 


U1 v2 U3 U4 U5 U6 U7 
Vi V2 V3 V4 U5 UE V7 
1 1 0 1 0 1 0] g 
H= 0 1 1 0 1 0 1 C2 
1 1 1 0 0 1 Oj g 
0 0 0 0 1 1 ija 
cy C2 C3 c4 


Figure 1. A sparse parity check matrix of Figure 2. A tanner graph of the left LDPC 
some LDPC code code 
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2.2. Genetic algorithm 

The GA is a type of evolutionary algorithm that belongs to the family of algorithms known as genetic 
algorithms. A GA’s population evolves via genetic operators inspired by biology’s evolutionary process [6], 
Darwin recognized that species evolution is driven by two processes: the process of selection and reproduction 
The reproduction of the fittest and most vigorous individuals is provided by selection, while reproduction is 
a phase in which evolution takes place. Travel salesman problem (TSP), job-shop scheduling problem (JSP), 
bandwidth-reduction problem (BRP), and linear ordering problem (LOP) are examples of permutation prob- 
lems. [7], which is a class of combinatorial optimization problems, the task is to arrange some genes 
(objects) in chromosome, with no duplicates, in a certain order that optimizes an objective function, where the 
representation of the chromosomes depend on types of the optimization problems [9], [10]. 

GA addresses the permutation issue by searching fast via the search space. It employs the selection, 
crossover, and mutation operators to produce superior chromosomes at the lowest possible cost [11]. The effi- 
ciency of using evolutionary algorithms to solve combinatorial optimization problems has been demonstrated 
(12}-[16]. It exists powerful algorithms, a nature-inspired algorithms like gaining-sharing knowledge based 
algorithm (GSK) [17]-[[19] which it has shown better results in solving optimization problems. The GA has 
several advantages such as: 

— Uses only the objective function’s evaluation, regardless of its nature (continuity, differentiability...), as a 
result of which there is more flexibility and a broader variety of applications. 


— Instead of a single iteration as in standard algorithms, generation adopts a parallel form by operating on 
several points at once. 
— Probabilistic transition rules (selection, crossover, and mutation probability) rather than deterministic ones. 

Many research have indicated that exhibiting a comprehension of the GA parameters’ interaction 
process, notably crossover probability, mutation probability, and population size, is the most important factor 
in evaluating the process. These factors are connected to each other in some way that impacts the GA efficiency. 
The optimal circumstance to use GA is when there is variety in the starting population with a high crossover 
chance and a low mutation probability [20]. 

It is important to note that the traditional crossover operator can not be applied to perform of per- 
mutation problems solution due to chromosomes arrangement of the genes is crucial, and no genes should be 
duplicated or missing [11]. Also, In comparison to other scenarios, it is more computationally expensive. The 
reason for this is that for offspring with duplicate numbers, a legalization step is necessary after each substring 
exchange. In such a case, the time required to complete a crossover operation increases fast as chromosome 
size increases, which can reduce the efficiency of permutation-based GAs [21]. Liu and Kroll in their research 
article developed a genetic algorithm did not use the crossover operator. It is important to note again, that 
GA has been used to find Automorphisms set for some block codes like bose—chaudhuri-hocquenghem (BCH) 
and quadratic residue (QR) codes of small length [23], also to compute the minimum distance of linear block 


codes [24]. 


3. GENETIC ALGORITHM-BASED METHOD 

In this section of the article, the genetic algorithm-based method is proposed, which uses an encoding 
that consists of treating an individual (permutation) as a sequence of numbers from 1 to the length of the code 
n. Also, these proposed method components work as explained in the next subsections. These components 
of the algorithm, which are the fitness function, which is used in the calculation of an individual’s fitness 
value, those fitness values are crucial in the choosing and construction of the individuals of the next generation 
through operators. The search space consists of n! individuals, each with n digits. The selection, crossover, 
and mutation operators will be explained and illustrated with figures. Then an overall organigram that shows 
how the algorithm works will be presented, identifying inputs and outputs. 


3.1. The search space and fitness function 

Let C be a binary linear code of length n, since our problem of finding the stabilizers set belongs to the 
optimization problems, the size of search space is linked to the code length n, This search space where the our 
proposed method will search, contains n! permutations. For all permutation o € Sn, each permutation will be 
associated to its corresponding permutation matrix, so every permutation of codewords will be in matrix form, 
including calculation of fitness values P,: 
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P= 0(In) (6) 
M,C C is a codewords set, such that, Vc; € Me: 


aiHT =0 (7) 


Ms, is matrix where its rows formed by codewords: 


Cy C11 C12 ess Cin 
C2 C21 C22 «++ Can 

Ms.=|.| =]. : : : (8) 
Ck Ckl Ck2 +++ Ckn 


applying the action of Sn on Ms,,.Vo € Sn s-t: 


Co(11) o(12) +++ Co(1n) 
Co(21) o(22) +--+ Ca(2n) 
a(Ms,) => Ms. Po = J . A P (9) 
Co(k1) Co(k2) +++ Co(kn) 
51 
$2 
o(Ms,)H* = | | | , where H is the sparse parity-check matrix (10) 
Sk 


The permutation of Mg, columns will generate another matrix of codewords if o(Ms,)HT = 0 (5). 
The selection of best permutations (individuals) will be based on the fitness values of permutations using the 
fitness function which is defined as follows: 


k 
fo = >_ wt(si) (11) 
i=l 


where s; is the syndrome of a codeword c; [25], and wt is the weight. The selection operator will need the 
values of each permutation which is calculated using the fitness function (11) in order to select that permutation 
or not. 

The Figure 3]shows the crossover operator which bases on the composition, which is chosen in order 
to ensure that all produced individuals within the search space and elements of S,, without relying on mutation 
due to the mutation operator probability of which is very low. Also, our method will use the mutation operator 
that consists a swapping of two gene’s position of an individual as figured in the Figure] this mutation type is 
chosen to enhance the convergence of the algorithm and to obtain new individual fitness of which are better. 


Pt = Parent 1 EsEAEVESEARSESES 
P2 = Parent 1 Ps titet7[2]{e6 [4] 5] 
3 2 4 


2 4 5 
ors [i ]e]s]]«]7]*] 5] 
ompiga-reo [a ]eTe]e]*]*]]7] 


Figure 3. Crossover operator 
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Chromosome EXEREE SEAS 
Mutated chromosome a HEEE GE 


Figure 4. Mutation operator 


3.2. The method inputs and outputs 


The following is how the GA-based method works: 
Inputs: 


— A codewords set Ms, 


The initial population size N; 


The number of generations Ng 
— The crossover probability pe 


The mutation probability p,, 


Outputs: 
— The Automorphisms permutations set 


The Figure Blis the genetic algorithm-based method organigram where the selection operator uses the 
fitness function values (6), and the stochastic crossover and the stochastic mutation operators are explained in 


Figures B]ana [4] 


Ms.» Pe: Pm, Ni and Ng 


Generating the initial 
population N=1 
Selection / Crossover / 
Mutation 
Fitness computation f, 


Store all o; S.t, So; = 0, 
then move on to next 
generation 


fo(Ms.) =0 
andN > N, 


Figure 5. Genetic algorithm-based method organigram 
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4. RESULTS AND DISCUSSION 

The results are obtained using parameters cited in Table The permutation is presented as a list 
where the positions are numerated from 1 to the length of LDPC code. Every error correcting code has an 
automorphisms group, therefore the set of automorphism permutations set exist for LDPC codes. Figure 6 
contains 160 automorphisms permutations produced by our GA-based method for [8,4,2] LDPC code and 12 
automorphisms permutations for [16,8,3] LDPC code listed in the Figure 7. 


Table 1. Paramters of GA-based method 


Parameter Value 

Initial population size 200 

Selection elitism 

Crossover rate 0.85 

Mutation rate 0.02 

Number of generations 30 
1,2,4,3,8,6,7,5] [1,2,4,8,3,6,7,5] [1,2,8,4,5,6,7,3] [1,2,8,5,4,6,7,3] [1,3,2,5,4,8,7,6] [1,3,5,6,2,8,7,4] [1,3,6,5,4,8,7,2 
1,4,3,2,6,5,7,8] [1,4,3,6,2,5,7,8] [1,4,6,8,3,5,7,2] [1,4,8,2,6,5,7,3] [1,5,2,8,3,4,7,6] [1,5,3,2,6,4,7,8] [1,5,3,6,2,4,7,8 
1,6,3,4,5,2,7,8] [1,6,8,4,5,2,7,3] [1,8,5,2,6,3,7,4] [1,8,5,6,2,3,7,4] [1,8,6,5,4,3,7,2] [2,1,3,4,5,7,6,8] [2,1,3,5,4,7,6,8 
2,1,4,8,3,7,6,5 2,1,5,3,8,7,6,4 2,3,1,5,4,8,6,7 2,3,7, 4,5, 8,6,1 2,4,1,3, 8,5, 6, 7 2, 4,7,3,8, 5,6,1 2,4,7,8, 3,5,6,1 
2,5,1,3,8,4,6,7] [2,5,1,8,3,4,6,7] [2,5,3,1,7,4,6,8] [2,7,3,5,4,1,6,8] [2,7,4,8,3,1,6,5] [2,7,5,3,8,1,6,4] [2,7,5,8,3,1,6,4 
2,8,1,5,4,3,6,7] [2,8,4,1,7,3,6,5] [2,8,5,1,7,3,6,4] [3,1,2,5,4,7,8,6] [3,1,4,2,6,7,8,5] [3,1,6,5,4,7,8,2] [3, 2,1, 4, 5,6,8,7 
3,2,1,5,4,6,8,7] [3,2,4,1,7,6,8,5] [3,2,4,7,1,6,8,5] [3,4,1,2,6,5,8,7] [3,4,1,6,2,5,8,7] [3,4,7,6,2,5,8,1] [3,5,6,1,7,4,8,2 
3,5,7,2,6,4,8,1] [3,6,1,4,5,2,8,7] [3,6,4,1,7,2,8,5] [3,6,7,4,5,2,8,1] [3,7,5,2,6,1,8,4] [4,1,2,8,3,7,5,6] [4, 1,3,2, 6,7,5,8 
4,1,3,6,2,7,5,8] [4,1,6,3,8,7,5,2] [4,1,8,6,2,7,5,3] [4,2,1,3,8,6,5,7] [4,2,1,8,3,6,5,7] [4,2,7,8,3,6,5,1] [4,2,8,7,1,6,5,3 
4,3,1,2,6,8,5,7] [4,6,1,3,8,2,5,7] [4,6,1,8,3,2,5,7] [4,6,8,1,7,2,5,3] [4,7,2,3,8,1,5,6] [4,7,8,2,6,1,5,3] [4,7,8,6,2,1,5,3 
4,8,7,2,6,3,5,1] [5,1,2,3,8,7,4,6] [5,1,2,8,3,7,4,6] [5,1,3,6,2,7,4,8] [5,2,1,8,3,6,4,7] [5,2,3,7,1,6,4,8] [5,2,8,1,7,6,4,3 
5,3,2,7,1,8,4,6] [5,3,6,1,7,8,4,2] [5,3,6,7,1,8,4,2] [5,6,1,8,3,2,4,7] [5,6,8,1,7,2,4,3] [5,7,2,3,8,1,4,6] [5,8,1,2,6,3,4,7 
5,8,1,6,2,3,4,7] [5,8,2,1,7,3,4,6] [5,8,2,7,1,3,4,6] [5,8,7,2,6,3,4,1] [5,8,7,6,2,3,4,1] [6,1,5,3,8,7,2,4] [6,1,5,8,3,7,2,4 
6,1,8,4,5,7,2,3 6,1,8,5,4,7,2,3 6,3,1,4,5,8,2,7 6,3,1,5,4,8,2,7 6,3, 7,5, 4,8,2,1 6,4,1,3,8,5,2,7 6,4,1,8,3,5,2,7 
6,4,3,7,1,5,2,8] [6,4,8,7,1,5,2,3] [6,5,7,8,3,4,2,1) [6,7,3,5,4,1,2,8] [6,7,5,8,3,1,2,4] [6,7,8,5,4,1,2,3] [6,8,1,4,5,3,2,7 
6,8,1,5,4,3,2,7] (6,8,4,7,1,3,2,5] [6,8,5,1,7,3,2,4] [6,8,5,7,1,3,2,4] [6,8,7,4,5,3,2,1] [7,2,3,4,5,6,1,8] [7,2,4,3,8,6,1,5 
7,2,4,8,3,6,1,5] [7,2,5,8,3,6,1,4] [7,3,2,4,5,8,1,6] [7,3,2,5,4,8,1,6] [7,3,4,6,2,8,1,5] [7,3,5,2,6,8,1,4] [7,4,2,3,8,5,1,6 
7,4,2,8,3,5,1,6] [7,4,3,2,6,5,1,8] [7,4,3,6,2,5,1,8] [7,4,6,3,8,5,1,2] [7,4,8,2,6,5,1,3] [7,5,3,2,6,4,1,8] [7,5,6,3,8,4,1,2 
7,5,6,8,3,4,1,2] [7,5,8,2,6,4,1,3] [7,6,3,5,4,2,1,8] [7,6,4,3,8,2,1,5] [7,6,4,8,3,2,1,5] [7,6,5,8,3,2,1,4] [7,6,8,5,4,2,1,3 
7,8,4,2,6,3,1,5] [7,8,5,2,6,3,1,4] [8,1,4,2,6,7,3,5] [8,1,6,5,4,7,3,2] [8,2,4,1,7,6,3,5] [8,2,5,7,1,6,3,4] [8,2,7,4,5,6,3,1 
8,4,1,2,6,5,3,7] [8,4,2,1,7,5,3,6] [8,4,6,1,7,5,3,2] [8,5,1,2,6,4,3,7] [8,5,2,1,7,4,3,6] [8,5,2,7,1,4,3,6] [8,5,6,1,7,4,3,2 
8,5,6,7,1,4,3,2] [8,5,7,6,2,4,3,1] [8,6,1,5,4,2,3,7] [8,6,4,1,7,2,3,5] [8,6,7,4,5,2,3,1] [8,6,7,5,4,2,3,1] [8,7,2,4,5,1,3,6 
8,7,2,5,4,1,3,6] [8,7,4,2,6,1,3,5) [8,7,4,6,2,1,3,5] [8,7,5,2,6,1,3,4] [8,7,6,4,5,1,3,2] [8,7,6,5,4,1,3,2 


Figure 6. Automorphisms set of [8,4,2] LDPC code 


[1, 2,3, 4,5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16] [1, 2,3,4,5,6,7,8,9,10,16,12,13,14,15,11] [1, 2, 3, 4, 5, 6, 7, 9, 8,10, 11,12, 13, 14, 15, 16] 
[1, 2, 3, 4, 5, 6, 7, 9, 8, 10, 16, 12, 13, 14, 15, 11] [1, 2, 3, 4, 8, 6, 7, 5, 9, 10, 11, 12, 13, 14, 15,16] [1, 2, 3, 4, 8, 6, 7, 5, 9, 10, 16, 12, 13, 14, 15, 11] 
[1, 2, 3, 4, 8, 6, 7, 9, 5, 10, 11, 12, 13, 14, 15, 16] [1, 2, 3, 4, 8, 6, 7, 9, 5, 10, 16, 12, 13, 14, 15,11] [1, 2, 3, 4, 9, 6, 7, 5, 8, 10, 11, 12, 13, 14, 15, 16] 
[1, 2, 3, 4, 9, 6, 7, 5, 8, 10, 16, 12, 13, 14, 15, 11] [1, 2, 3, 4, 9, 6, 7, 8, 5, 10, 11, 12, 13, 14, 15, 16] [1, 2, 3, 4, 9, 6,7, 8, 5, 10, 16, 12, 13, 14, 15, 11] 


Figure 7. Automorphisms set of [16,8,4] LDPC code 


To be mentioned, each combination of two automorphisms permutations is an automorphism permu- 
tation, if the set contains all generators of automorphisms group, then we can obtain the others automorphisms 
permutations easily. The Table[2]shows statistical measures of 32 runs of GA-based method for [8,4,2] LDPC 
code, which shows the efficiency of our method for finding an important automorphisms set, in some runs, we 
get an important number of automorphisms in few number generations (set of 160 Automorphisms permuta- 
tions in 8 generations). 


Table 2. The statistical measures 
Mean Median Standard deviation Best Worst 
151.09 152.5 11.09 160 104 


Indonesian J Elec Eng & Comp Sci, Vol. 25, No. 2, February 2022: 1059-1066 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 o 1065 


5. 


CONCLUSION 
In this paper, the genetic algorithm-based method has been proposed for finding an important auto- 


morphisms set of a given LDPC code, which can be used in improving their decoding algorithms (the hard 
decision algorithm and the soft decision algorithm). It showed good results for LDPC codes in the short block 
length regime. Our future work is to optimize our method and its genetic parameters, and combining it with 
the GSK algorithm in order to process long LDPC codes. 
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